NIKHEF-H/92-18 



The renormalization of the axial anomaly 
in dimensional regularization 



S.A. LarinQ 

NIKHEF-H 
P.O. Box 41882 
1009 DB Amsterdam. 



Abstract 

The prescription for the 75-matrix within dimensional regularization in 
multiloop calculations is elaborated. The three-loop anomalous dimension 
of the singlet axial current is calculated. 
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Since the invention of dimensional regularization [Q] and the minimal sub- 
traction (MS) scheme [0] a lot of attention was paid to the problem of the 
75-matrix within dimensional regularization. The following approaches to this 
problem were used in practical calculations: the prescriptions based on the orig- 
inal definition by 't Hooft and Veltman ]!]] || |^], keeping the four-dimensional 
anticommutation relation for 75 in D-dimensions [|| and dimensional reduction 
Q. Discussions of the 75-prescriptions can be found e.g. in (Tj ||- 

The effective approach to perform multiloop calculations involving the non- 
singlet axial current in dimensional regularization was developed in |]] [10] 



where deep inelastic sum rules were calculated up to (and including) the three- 
loop level in QCD. The effectiveness of this approach is confirmed by its use in 
such an advanced calculation as the calculation of the deep inelastic structure 



function F3 at the two- loop QCD level [11] 



In the present paper this approach is elaborated for the cases of the pseu- 
doscalar current and the singlet axial current. The three- loop anomalous di- 
mension of the singlet axial current is calculated by imposing the requirement 



that the axial anomaly relation [ITI3] [13] should preserve the one-loop character 



[14] in dimensional regularization. 

Throughout the paper we use the MS'-scheme (2| or its standard modifica- 
tion, MS'-scheme jnj, to perform renormalizations. The dimension of space- 
time is defined in the standard way as D = 4— 2e. All calculations are performed 
within massless perturbative QCD. 

1. The non-singlet axial current. Let us first consider the case of the 
non-singlet axial current: 

Jl a (x) =^(x) llll5 t a ij(x), (1) 

where ip is a quark field and t a is a generator of a flavor group. 

In our opinion, the most practical definition of 75 for multiloop calculations 
in dimensional regularization (and the only one known to be self-consistent) is 
the original definition due to 't Hooft and Veltman Q: 

75 = *^J £ t'li;2t'3^47^l7z'27z'3T / V4 ) (2) 

here the Levi-Civita e-tensor is unavoidably a four-dimensional object and 
should be taken outside the /^-operation where any object can be safely con- 
sidered as a four-dimensional one; the indices v\ . . . 1/4 are D-dimensional inside 
the i?-operation as all other indices within dimensional regularization. But 75 
defined by eq.(H) does not anticommute anymore with the L>-dimensional 7^. 
That is why in order to define the axial current correctly one should use (see 
below comments after eq.(|lC)D) the symmetrical form of the axial current: 

= 2^(7^75 - 757m)* < V' 75 = ijfv 1 v 2 u 3 v i lv 1 lv 2 lv^v A - (3) 

In principle it is possible to perform the calculations using this definition of the 
axial current. But one can simplify the definition drastically. Let us commute 
7^ in the first term in (^) to the right. The D-dimensional metric tensors g tiVi 
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arising during commutations can always be taken outside the ii-operation where 
they can be safely contracted with the e-tensor as four-dimensional objects. So 
we receive the following definition of the non-singlet axial current: 



This is exactly the prescription proposed in ||. Thus we proved equivalence of 
the definitions (||) and (||) within dimensional regularization. To be sure that 
no holes are missed in this general proof we computed the axial vertex (see 
below eq.(||)) at the three-loop level, using both definitions. The results are 
identical. But the definition (Q) is more compact and saves computational time 
enormously so we will use it in what follows. 

Since the anticommutativity of 75 is violated by definition (||) , the standard 
properties of the axial current and Ward identities which are valid e.g. in such 
a basic regularization as the Pauli-Villars regularization are also violated. In 
particular the renormalization constant Z^f s of the non-singlet axial current in 
the M/S-scheme is not equal to one any more. It was calculated in the three-loop 



approximation in QCD in [10]; we remind here the two- loop expression: 

1 22 4 

Ztl s = l + a 2 - e {-C F C A --C F n f ) : (5) 
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where we use the notation a = ^ = j-Jr-y for the strong coupling constant, Cf 
and Ca are the Casimir operators of the defining and the adjoint representations 
of the color group and rif is the number of quarks flavors. The relation between 
renormalized and bare operators is (O)r = Z (0)b- 

To restore the renormalization invariance of the non-singlet axial current 



(i.e. to nullify its anomalous dimension) one should perform [16] an extra finite 
renormalization or in other words to introduce the extra finite renormalization 
constant Z£ s (a). So the correct renormalized non-singlet axial current is: 

1 

3? 



(4 a ) R = Zr(a) Zfi s (a) (fp B = Z% s i^e^^^^t^B, (6) 



where ip B = Z^ is a bare quark field. The full anomalous dimension can be 
now nullified: 



if (a) = log(Z 5 - Zlfs) = 0(a)— + log = 0. (7) 



Using this equation one can obtain Z£ s from the given ZV^a. But since the 
renormalization group /3-function starts with an a 2 -term one can obtain Z£ s (a) 
only in the approximation which is one order in a less than the given approxi- 
mation of Zj^ s . That is why it is better to use the recipe of |J and to find Z£ s 
from the requirement that the renormalized axial and vector vertices coincide: 

Z™R MS < i> 4 a (0) V >= Rms < ? ^(0) V > 75, (8) 

where J^(x) = ip(x)j fJi 'f5t a "ip(x) is the vector current. This relation means that 
anticommutativity of the 75-matrix is effectively restored, so the standard Ward 



3 



identities are also restored. This prescription for Z§ s automatically ensures zero 
anomalous dimension (0) because the anomalous dimension of the vector vertex 
is naturally zero. 

The first impression is that the calculation of the axial vertex is rather 
cumbersome within adopted prescription (Q) for the axial current: taking the 
e-tensor outside -R-operation seems to leave inside R-operation three extra un- 
contr acted indices ^1,^2,^3. But it is always possible to contract first a quan- 
tity under consideration with an extra e-tensor to produce a scalar object. In 
our case we first multiply the axial vertex with 7^75 = p 1 p 1 p2p- i l pil P2I p-j, arm 
take the trace. Now outside the R^s-operation the product of four-dimensional 
e-tensors can be represented in a standard way as the determinant of the four- 
dimensional metric tensors: 

£pu 1 l/ 2 ^3 £ fJ-PlP2P3 = 9v 1 p 1 {9l/2P29lS3p3 ~ 9v2P2,9v 3 p2) + ••• (9) 

Then outside the -R^zs-operation these four-dimensional metric tensors can be 
safely considered as D-dimensional ones (which will add only inessential 0(e)- 
terms to the renormalized axial vertex). These D-dimensional metric tensors 
9/iipj can be safely taken inside the Reoperation. So one obtains finally a scalar 
expression inside the Reoperation containing only L>-dimensional objects which 
makes the practical calculations straightforward. 

Calculating the axial vertex and the vector vertex in eq.(|8|) one finds 
The three- loop approximation for Z£ s was obtained in [10]; we recall here the 
two-loop expression: 

107 2 

Zl s = 1 + a(-4C F ) + a 2 (22C 2 F - —C F C A + -C F n f ). (10) 

y y 

The independence of Z£ s on the log(p 2 / /j?) (where p is the momentum of quark 
legs in the vertices in eq.(||)) gives a strong check of the whole prescription. For 
example, trying to use the axial current in the naive form ([!]) with 75 defined 
in one would obtain that log(p 2 /// 2 ) does not cancel in Zg s which excludes 
the possibility to use this naive form. 

Note that both Z^f s and Z^ s are gauge independent quantities (which was 
checked by calculations in an arbitrary covariant gauge). We would like to note 
also that the finite constant Z£ s (a), like the usual ultraviolet renormalization 
constants, does not depend on the choice of the concrete modification of the 
MS- like schemes: whether it is calculated within the MS'-scheme itself or MS- 
scheme or G-scheme [p"7|]. 

So to calculate any quantity involving the non-singlet axial current one can 
use the prescription (||) with Z^s an d Z$ s given in (||) and (|Toj) . This pre- 
scription has all typical features, see @, of the approach developed in Q, where 
D-dimensional indices split into 4-dimensional and {D — 4)-dimensional indices. 
But in our approach all indices during the calculations are D-dimensional which 
avoids in practical calculations all complications connected with splitting in- 
dices. 

In the case of several 75-matrices in one fermion line one can naively an- 
ticommute them (because of the validity of eq.rtS)) as in the prescription with 
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anticommuting 75 M and use the standard property 75 = 1. So e.g. the corre- 
lator of two axial non-singlet currents automatically coincides in the approach 
under consideration with the correlator of two vector non-singlet currents. 

2. The pseudoscalar current. Let us apply now this prescription to the 
case of the pseudoscalar current: 

P(x) = ip(x)j 5 tp(x) = ^£^2^4*07*17^7^7^ (11) 

We will not distinguish singlet and non-singlet cases for the pseudoscalar current 
because the pseudoscalar current does not generate closed fermion loops in 
the massless case. As in the case of the non-singlet axial current, to obtain 
the correct renormalized pseudoscalar current we should introduce the finite 
constant Z P in addition to the usual ultraviolet renormalization constant Z P MS 
in the MS-scheme: 

(P) R = Z P (a)Z p MS (a)(P) B = Z%Zl Is i-£ UlU2U3U4 ^ B i Vl j U2 j U3 'y 1/4 tl)B (12) 

We can calculate first in the standard way the renormalization constant Z P MS 
within the MS-scheme. In all calculations we use the program Mincer ||i~8fl 
written for the symbolic manipulation system Form [19]. This program com- 



putes analytically the three-loop massless diagrams of propagator type which 
is sufficient to calculate any renormalization constant within MS-scheme at the 
three-loop level. 

The three-loop approximation for Z P MS in the MS-scheme is: 

Z P ms = ^ + ^-C F l)+a 2 [C F C A (^ + ^ + C F n f (-^-f e ) + CH^- 2 - 

+n*\r r w 44 + 110 4- 8C3 58 urr^ 121 257 599 ) 

+a ^^9? + ^ + — -^e ) + CfCa{ ~9^ ~ 54? - W 



2 4 22 17 2 33 215 3203, 

^ ( ~9? " 27? + 27? + CfCa ^ ~T2? + ^6T' 



where (3 is the Riemann zeta-function (£3 = 1.202056903...). To restore the 
Ward identities we can, as in the case of the non-singlet axial current, define 
the finite renormalization constant Z P from the requirement of coincidence of 
the pseudoscalar and scalar vertices: 

Z P R MS < ? P(0) V >= Rms < ^ V^(0) V > 75, (14) 

So the anticommutativity of the 75-matrix is effectively restored and it is an- 
ticommuted out of the pseudoscalar vertex. Calculating the three-loop pseu- 
doscalar and scalar vertices we find the three- loop approximation for Z P : 

2 4 64 856 

Zl = l + a (-8C F ) + a\-C F C A + -C F n f )+a 3 [C F C A n f (—C 3 + —) 



+( H 1-208(3 ~ 9 -§)+ l ^C F n) + C F C A (608{ 3 - ^ 
64^ 580, „o, _ 304 



+C F n f (- Y ( 3 - —) + C^(-384C 3 + — )] (15) 



Again the cancellation of log(p 2 //i 2 ) in Z P provides a good check of the cal- 
culations. Another check is the gauge independence of both Z P MS and Z P (all 
calculations were done in an arbitrary covariant gauge). Note that Z P for the 
pseudoscalar current differs from Z§ s for the non-singlet axial current. 

Note also that the full renormalization constant Z P Z P MS of the pseudoscalar 
current does not coincide with the renormalization constant of the scalar current 
Z^ but their anomalous dimensions do coincide: 

f A log(Zf Z P MS ) = ^ A log Z ^ = +a(3C F ) + a 2 (\c 2 F + ^C F C A - \c F n f ) 

,„ 3[ 129 3 129 2 ,11413 2 

278 35 

+C|n/(24C3 - 23) + OC^/^Ca - ^) - —C F nJ]. (16) 

We would like to note that this our result agrees with the original calculation 
of the three-loop anomalous dimension of the quark mass in the MS'-scheme in 
pQl and provides thus the independent check of that calculation. 

Thus one can use in all calculations involving the pseudoscalar current the 
prescription fll2| ) with Z P and Z P MS given in ( fl5| ) and (ff^) at the three-loop 
level. In the case of several pseudoscalar vertices in one fermion line one can 
(since the anticommutativity of the 75 is effectively restored by the prescription 
(0)) naively anticommute 75 and use the standard property 75 = 1. So e.g. the 
correlator of two pseudoscalar currents automatically coincides in the considered 
approach with the correlator of two scalar currents. 

3. The singlet axial current. Let us consider now the case of the singlet 
axial current which we define in the analogy with the non-singlet current (||): 

Jl = V^TsV' = i^w^v^lvxlvfiv^- (17) 

It is known that the singlet axial current is nontrivially renormalized be- 
cause of the axial anomaly and the renormalization constant of the singlet axial 



current is nontrivial at the two- loop level [ 12 1 , pi] ] . To receive the correct renor- 
malized singlet axial current we need, as in the previous cases, to introduce 
the finite renormalization constant Z§ in addition to the standard ultraviolet 
renormalization constant Z S MS within the MS-scheme: 

(^)-R = Z 5 Z ms( j ^)b = ZlZ s MS i—e^ UlU2U3 i> B j Ul -/ U2 -/ U3 ip B - (18) 

We can calculate within the MS-scheme the renormalization constant Z S MS of 
the singlet current at the three-loop level: 

Z' MS = 1 + a 2 [C F C A (f e ) + C F n f (^-)) + a 3 [C F C A n f (-^ + ^) 

_ 484 3578. „ 2 . 20 26. ^ 2 , 308. ^ 2 . 22.. 

+CfC* a (-— 2 + —) + C F n%— 2 + — ) + C\ C A (-— ) + C F n f (--)Hl9) 

Now the problem is how to fix the finite renormalization constant Z\ . As in the 
previous cases one should restore within dimensional regularization the stan- 
dard properties of the singlet current which exist in such a basic regularization 
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procedure as the Pauli-Villars regularization. We can not anymore impose for 
this purpose the coincidence of the axial and vector vertices (||) because the 
singlet current generates closed fermion loops and we cannot anticommute 75 
outside the singlet axial vertex. 

To fix Z§ one can require within dimensional regularization the conservation 



of the one- loop character [14] of the operator relation of the axial anomaly which 



is valid in the Pauli-Villars regularization: 

(d,4) R = a^(GG) R , (20) 

where GG = e^G^G^ and G% = d^A* - d v A^ + gf abc A b ^ is the gluonic 
field strength tensor. 

Let us consider in detail renormalizations of both sides of the anomaly 
relation. The divergence d^J^ is renormalized multiplicatively in the same way 
(|lg|) as the current itself: 

(d.J^n = ZIZ s MS (d^) B . (21) 

But the operator GG mixes under renormalization: 

(GG) R = Z GC ,(GG) B + Z GJ (d^) B . (22) 



It is known [22] j23[ that to explain why d^A^ does not mix under renormaliza- 
tion and GG does mix it is convenient to represent GG as the divergency of the 
axial gluon current: 



GG = 



— ^ £ ^uiv 2 u 3 {A^ i d U2 A'l i + -gf^A^A^A^). (23) 

The current Kn is not gauge- invariant. That is why the gauge invariant current 
in some "good" gauges (e.g. axial gauge or background- field gauge) can not 
mix with the gauge- variant operator K^. Then the divergences d^J^ and GG of 
the current and are known to be renormalized as the current themselves. 
So d^Ap does not mix with GG in these "good" gauges either. But since both 
divergences are gauge invariant this non-mixing is valid in any gauge. 

To understand the restrictions on the renormalization constants Z G q and 
Zqj it is useful to take the renormalization group divergence of eq . (p2^) : 

li 2 -^{GG) R = l GG {GG) R + 7gj(VJ)*, (24) 

where anomalous dimensions are defined for the case of operator mixing as 
follows: 

d dz {1) 00 z (n) (a) 

(Oi) R = ZijiO^B, 7ij = (n^ZikXZ- 1 )^ = -a-£-, Z i5 = 1 + E ^T^( 25 ) 

" n=l 

From the renormalization invariance of the anomaly: 

M 2 ^(VS)* = M a ^a^(G6)«, (26) 
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one receives now: 

Ij^JDr = a^fd + i G a){GG) R + 7 gj(^J>] = + 7 GG + ^7gj)(^J>. (27) 

From this equation one can assume the restrictions on the anomalous dimen- 
sions: 

7 G g = -^, 7GJ = («f )-S5- (28) 

Strictly speaking eq.(|27|) itself permits a more general solution than the solution 
given in eq.(|2^). But the direct calculation at the two- loop level (see below 
eq . (|3l|) ) supports eq.(p8|). 

To calculate Z| one can calculate matrix elements of the l.h.s. and r.h.s. 
of the anomaly operator equation (p0|) between gluon states (so the famous 
anomalous triangle one-loop diagrams plus higher loops appear in the l.h.s.). 
To be more precise we calculate for the l.h.s of eq.(^) the following Green- 
function: 

< A d^Jl A >= 

= RJTS / d'xd'ye^y < T^(*)d M J»^(0) "",(29) 

where 'amputated' means that one-particle-irreducible diagrams with ampu- 
tated external gluon legs are considered. As it was explained above, the essential 
point of the calculation is that the product of two e-tensors can be substituted 
by the determinant of metric tensors which can be taken as D-dimensional ones 
inside -R-operation. Some typical three- loop diagrams contributing to eq.(p9|) 
are shown in fig.l. 

The result of the three-loop calculation in the MS-scheme is: 

1 323 
< AdptfA >= 24n f a{l + a[C F (4) + C A (6 + 2£ - -f)] + a 2 [C F C A (— + 8£ - 2f) 

+CFn f (~ + 8C 3 ) + C|(-6) + C A n f (-^ - 12 C 3 + f £ - 5 -f) 
^2/4537 OA „ 1„. 2 2467 A 13K 2 7^ 3 Ms „ , nn . 

+c\(— + 4 C 3 + 3 Cs e - -^e + ^ - -^e - - 8 e + T6 e)Mm 

where £ is the gauge parameter in an arbitrary covariant gauge so the gluon 
propagator is (g^ — £,^r-)/q 2 ■ We omitted in this result terms with \og{p 2 / n 2 ). 



For the r.h.s of the anomaly equation (20) we calculate at the two- loop 
level the analogous matrix element, as for the l.h.s (]29|). Some typical two- 
loop diagrams contributing to the r.h.s can be obtained from the diagrams of 
fig.l by shrinking the upper fermion loop into a point. Renormalization is 
done according to (|22|). The necessary approximations for the renormalization 
constants are: 

z gg = l + <4(-yCU + !"/)] 

2r l. 44 121 2 4 2 , K 5^ 17^ 2 M 

+a [^(— g° An f + ~Y C A + g n /) + -(°F n f + 3 C ^ n / - Y° a)] ' 

Z G j = a-\2C F . (31) 
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The validity of the restrictions of eq.(p8|) is confirmed in this approximation. 
The two-loop result in the MS'-scheme is: 

a _L <A GGA> = 24n f a{l + a[C A (6 + 2Z--f)]+a 2 [+C F n f (-— + 8C3) 

+C A n f (-^ - 12 C 3 + f e - + C£(^ + 4Cs + 3C 3 € 
1 2 2467 a 131 2 7 A o 3 , , 

One can see that without the finite renormalization constant Z| the l.h.s of the 
anomaly relation (|30|) and the r.h.s. (^) do not agree within the MS-scheme. 
We can obtain the desired finite constant Z| restoring the one-loop character 
of the anomaly relation, i.e. dividing ( |3~2"1) by (|30|): 

1 D7 31 

Z s 5 = 1 + a(-4C F ) + a 2 (22C 2 - — C A C F + — C F n/). (33) 

y is 

We should stress that the difference between the singlet constant Z| and the 
non-singlet constant Z£ s is only in the C F nf -term. This difference is due to 
only the light-by-light-scattering type diagrams (the type shown first on fig.l). 
The independence of the obtained finite constant Z| on log(p 2 /// 2 ) gives strong 
confirmation that both sides of the anomaly relation are really matched. 

Thus both sides of the axial anomaly relation receive the non-trivial higher 
order corrections if one considers their matrix elements. But these corrections 
are matched and the one-loop character of the operator equation (which is valid 
in the Pauli Villars regularization) can be preserved in dimensional regulariza- 
tion. 

Now we can calculate the full anomalous dimension of the singlet current in 
0(a 3 ) approximation: 



75(a) = /i 



2 d)og(Z§Z' MS ) 



+a 2 (-6C F n f ) + a 3 (-^C F C A n f + \c F n) + 18C F n f ). (34) 



The first term agrees with the calculation [21]. 

It is interesting to consider now the transition to the QED case in the 
result ([k]) for the matrix element of the divergency d^J^ between gluon states, 
multiplied finally by the finite constant Z|. The transition from the QCD 
case to the QED case can be done by simple substitutions: Ca = 0, C F = 
1, nj/2 = rif, a s = a. Making this substitutions we find that the only surviving 
contributions to the matrix element of d^J^ between two photons are the famous 
one-loop triangle diagram and the three-loop diagrams of the light-by-light- 
scattering type (the first type on fig.l). The fact that these light-by-light- 
scattering diagrams give a non-zero contribution to the matrix element of <9 M 
between two photon states was discovered originally by direct calculation in 



1 24]. These three- loop diagrams give the correction to the width of the neutral 



pion decay into two photons. 
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Thus in all calculations involving the singlet axial current one can apply the 
calculational power of dimensional regularization by using for the singlet ax- 
ial current the prescription (Jl8|) with corresponding renormalization constants 



given in (|l9|) and (j33|) 



The generalization of the considered 75-prescription for the massive case is 
straightforward. Within the MS-scheme the ultraviolet renormalization con- 



stants do not depend on masses [25|. The same is valid for the finite renormal- 
ization constants Zr y . so the obtained finite constants can be applied also in the 
massive case. 

I am grateful to the collaborators of the Theory group of NIKHEF-H for helpful 
discussions. 
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